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Abstract 

We consider the "thermodynamic limit" of a d-dimensional lattice of hyperbolic dynamical systems on the 2-torus, interacting 
via weak and nearest neighbor coupling. We prove that the SRB measure is analytic in the strength of the coupling. The proof 
is based on symbolic dynamics techniques that allow us to map the SRB measure into a Gibbs measure for a spin system on a 
(d + 1) -dimensional lattice. This Gibbs measure can be studied by an extension (decimation) of the usual "cluster expansion" 
techniques. 
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1 Introduction and Main Results 



In recent years a lot of attention has been devoted to the relation between non equilibrium statistical mechanics and 

dynamical systems theory. According to the point of view of Ruelle, Cohen and Gallavotti [Ru78] [GC95] a mechanical 
system evolving in a steady state can be described by a hyperbolic dynamical system and its properties can be deduced 
from the "natural" or SRB distribution (see below for a precise definition) associated with this dynamical system. This 
line of investigation has already produced several interesting results both analytical, like the "Fluctuation Theorem" , 
see [GC95], or numerical, like the works of Evans and Morris, see [EM90], and Moran and Hoover, see [MH87]. 
Nonetheless, most of the work has been devoted to low dimensional dynamical system, due to their accessibility both 
to analytical and to numerical study. In this paper we want to study the properties of the SRB distribution for a class 
of simple systems in very high dimension. For more references on this kind of systems see [JP98] . The precise model 
model we study here is taken taken from [BKL03]. 

Wc start considering a linear hyperbolic automorphism of the 2-torus T^. To be definite, we will always consider 
the so called Arnold cat map Sq -.T"^ defined by the action modulus 27r of the matrix 

i)^ (1) 

Note that the matrix A admits two orthogonal eigenvectors v± whose respective eigenvalues A± are such that A+ > 
1 > A_ and A+A_ = 1. For this reason the dynamical system sq is uniformly hyperbolic and the stable and unstable 
manifolds at any point <^ G are given by W^{t) = (f) + v±t mod27r. 

^From So we can construct the uncoupled lattice dynamics by considering as phase space the Cartesian product 
T = (T )^ (namely any point ip ^ T has Cartesian components {ip^j^^j^d), equipped with the metric d{ip,'4'') = 

2~I^I^I(^/)J, V'^) where d{ijj^,tfj'^) is the usual metric on and |^| = ^r ^ G Z**. On T the map So acts 

simply as 

5«(^) = {SoW)^ = so(V?). (2) 

Note that the stable and unstable manifold of Sq at a point -0 arc the Cartesian product of the stable and unstable 
manifold of sq for the points G T^, i.e. W^^{Q = ip + J2(^ '^ofiCc T^od 27r, where Wofi is the tangent vector to T 

that has null component on the tangent space to every but for where it coincides with v± . The action of Sq on 
W^{'tp) is naturally given by a diagonal linear transformation. 

We observe that the special choice of the matrix A plays no role in the following. Indeed we will show in Appendix A 
that our results stay true if we replace sq with any uniformly hyperbolic analytic automorphism of T , not necessarily 
linear. 

To add a coupling to this system we consider an analytic function g -.T and define 

5|(V) =^so(Vi) + £5(P^V')'= so(V«) + sfHi') (3) 

where (p^ip)v ~ V'r;+j, *-e. p is the group of the translations on ^Z". This means that the function f : T ^ T, whose 
^ component is — g o pi ^ is translation invariant. We want / to be short ranged: let the nearest neighbor sites of 
the site ^ be nn{^) = {rj : |^ — ?7| < 1}, we will assume that g depended only on ipnn^o)^ where wc have used the 

notation tjjv = {V-'cIC G ^} foi' ^ C 'Zf' . This implies that S^(^JJ) depends only on ijjnn(^)- More generally we could 
have assumed that g depends only on ipv where V is any finite subset of if containing but this would not have 
changed the substance of the following arguments. Moreover we will take g analytic in all its arguments. 

The dynamical system Sg admits many invariant measures. Among them is the "natural" or SRB measure defined 
as the weak limit of the volume measure on T under the evolution defined by S^, when such a limit exists and is ergodic. 
Being T infinite dimensional, to properly define this concept we will consider finite dimensional approximations. Let 
T/v = (T^)^™ where Va^ is the cube of side 2 A'' + 1 in if centered at the origin. To define the restriction of to 
T/v we have to fix the boundary conditions: we choose periodic ones. To this extent note that T/v can be naturally 
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identified with the submanifold of T formed by the points periodic of period 2N + 1. Moreover leaves such a 

def I 

manifold invariant so that we can define S^^n = ^s\t-h' confusion can arise we will suppress the index A'^. 

We can now define the SRB measure for S^.n as ^ifj^^ = limr^oo T'Yiit=Q i^s,Ny lA^ where the limit must be 
understood as a weak limit and yP^ = HjeVN '^^i/i'^^)'^ ^he Lebesgue measure on Tn. The existence of such 
a measure follows from rather general theorem on hyperbolic dynamical systems, if e is sufficiently small (sec for 
example [GBG03] and references there). Moreover jif/^^ is ergodic, always for small e, and we have that ^fp^{0) = 

IJ,% almost every ^, where O is an observable, i.e. a Holder continuous function from 
Tjv to ]R. This means that fJiff'^ is the statistic of S^- It is well known that the SRB measure is still well defined in the 
limit N —>■ oo, for e small enough. This was first proved by Bunimovich and Sinai in [BS88]. Starting from this work, 
the model eq.(3) (or similar models of coupled expanding automorphisms of the circle) has been widely studied in the 
literature, see for instance [PS91],[BK95],[BK96],[BK97] [JM95], [JP98], [BEIJK98], [JdLOO]. Many properties of such 
systems are well known, mainly uniqueness of the SRB state in the thermodynamic limit and exponential decay of 
correlations (see [BK95],[BK96],[BK97] for a proof of these properties). 

We further investigate the regularity properties of the limiting measure. Wc show that /i^^^ depends analytically 
on e. This means that if we consider an analytic observable O, i.e. an analytic function from TJv to M, we have that 
nfP'^ [O) is an analytic function in a domain that depends on the analyticity properties of O. 

The main point of this work is to show that such a property remains true when — » oo, i.e. we want to show that 
the domain of analyticity of ii^^ does not shrink to when ^ oo. More precisely we say that O : T ^ R is a 
local observable if it depends only on ipv for some finite y C Z . We can summarize our main results as follows 

Theorem: Given as above and a local observable O we have 

(1) fjb^^^{0) = lim;v->oo Mw^^(^) exists uniformly in s for e < eq independent on O. 

(2) if O is local and analytic then ij.^^^{0) is analytic in s for e < eo{0). 

The proof is mainly based on the possibility of mapping the SRB distribution into the Gibbs state of a suitable spin 

system on Z ^ and on the extension of classical techniques used to study such Gibbs states (i.e. cluster expansion) 
to the particular ones that occur in our system. The key point in order to get analyticity of the measure is proving 
that the SRB potentials {i.e. the potentials of the Gibbs state the SRB measure is mapped into) are rapidly decaying. 
Once this decay is proved, analyticity follows via standard techniques. Analyticity of the measure and convergence of 
cluster expansion imply in particular uniqueness of the Gibbs measure and exponential decay in space and time of the 
correlations of Holder continuous observables, see for instance [GBG03]. Our proof can also be adapted to the case of 
coupled analytic expanding circle map: in fact also these models can be mapped into spin systems, and proceeding as 
below one can prove that the SRB potentials satisfy the same decaying properties. 

The rest of the paper paper is organized as follows. In section 2 we give a brief review of the main properties of 
smooth uniformly hyperbolic systems and we briefly describe the construction that allows the above quoted mapping. 
The detailed proofs of this properties are postponed to sections 3, 4 and 5. Finally in section 6 we complete the 
proof of the main theorem. Appendix A contain a direct extension of our results to the case in which the uncoupled 
dynamics is not linear. Appendix B deals with an application. In the contest of the physical application of dynamical 
systems (see the beginning of this introduction) a special status has been given to a particular observable, the phase 
space contraction rate defined as ri+{ip) = log | det(£'5£('^))| where DSg is the differential of S^. Being our system 
infinite it is more interesting to study the local phase space contraction rate rjv {ijj) defined by taking the determinant 
of a (large) minor of DSg. We show, for a large class of couplings /, that rjv has a positive average and that it obeys 
a large deviation principle, i.e. its large deviation are asymptotically described by a free energy functional. 
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2 Anosov Systems 



2.1 Geometric Properties 

A dynamical system on a smooth compact manifold, whose dynamics is given by a uniformly hyperbolic invertible 
map is called an Anosov system,. /,From the general theory wc know that Anosov systems arc structurally stable, 
namely, given two Anosov diffeomorphisms S, S' on a manifold Q that are sufficiently close in the topology, there 
exist a conjugation H : fl< — >0 such that S o H = H o S'. 

In our situation this implies the existence of a map : 7jv< — such that 

Seoh, = h,oSo, (1) 

at least if s is small enough (a priori not uniformly in N). The first step of our proof consist in showing that is 
analytic in e uniformly in N. More precisely we will construct directly for N = oo. Its finite N version can be 
obtained by restricting it to Tjv- We note that is, in general, only Holder continuous in the variable ip. By this we 
mean that there exist constant c and (3 such that d{hE{tp),hs{ip')) < cdi^ipjtp')^ . For this reason we cannot say that 
the SRB measure of is just the image under the map hg of the SRB measure for Sq, i.e. of the Lebesgue measure 



on T^. Notwithstanding will play a crucial role in the construction on ^fj 

As we saw in the introduction the tangent space T^T to T on a point -0 can be split in two subspaces and E'^ 
such that T^T = E'^ © E^. Moreover the distributions E^ are continuous and invariant under Sq, i.e. [DSq = 
Eg^^ and we have 

\\DSSw\\<CXl\\w\\ iovw&E^ 
||£>5o-"w|| <CA;"||u;|| for u; g £+ (2) 

E^ and E^ are called the stable and unstable subspaces, respectively. In the case of 5*0 all these properties are trivially 
true. In particular we can consider on E^ the basis formed by the vectors {wof±}jg^<i- 

We will show in section 4 that such a splitting can be constructed also for Se, again uniformly in A'', i.e. we will 
prove the existence of the stable and unstable subspaces Ef,^ for 5*^. Moreover we will show that is an analytic 

function of e, although is only Holder continuous in if). This will turn out to be the right regularity to study the SRB 
measure. To do this we will directly construct the vectors of the basis coincide with {wo^j_}^g^d 

for £ = 0. 



SRB 



2.2 Symbolic Dynamics 

The main property that allow to study analytically the SRB measure for an Anosov map S acting on a manifold 
M. is the existence of Markov partitions. We call a collection Qj, i = 1, . . . , n, of closed subsets oi M. & partition if 
UiQi = ^A and Qi n Qj = dQi n dQj for every i ^ j. For every sequence a = {<Jt}teZ ^ {1j • • • ; f^}^ can define the 
set X{a) = n^_(^5'*(Qcrt)- Due to the hyperbolicity properties of S if Qi are small enough X{a) contains at most 
one point. This allow us to construct a symbolic dynamics i.e. a map from a subset E of {1, ... , n}^ to A^. In general 
the structure of the subset E is very complex but for Anosov systems it is possible to construct particular partitions 
for which the set S can be described easily. Given a partition Q we call compatibility matrix the nxn matrix C given 
by Cij = 1 if int(5Qi) n int(Qj) ^ and otherwise. We say that Q is a Markov partition if the set S is formed by 
the sequences a such that Cai,ai+i = 1 for every « G Z. This means that the sequences that satisfy the above nearest 
neighbor condition code all the points of A^. In such a case we will denote £ = {!,..., n}§. 

We now show how to construct a Markov partition for our model. We start with sq. A Markov partition Q = 
{Qi, i = 1, . ■ ■ ,n} for So acting on T can be easily constructed starting from its stable and unstable manifolds. Such a 
construction is standard and can be found, e.g. , in [GBG03]. Let C be its compatibility matrix and cq the associated 
symbolic dynamics. 
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It is important to note that cq is Holder continuous in the sense that there; exist constants c and [3 such that, 
for any two sequences a, a' G {!,..., n}p, d{co{a),co{a')) < cd{a, cr'Y , with d[a, a') = e"^^'^''^ ^ where #(cr, a') the 
biggest integer such that Oj = a'j, V|j| < #((t, a'). In this case we can take /3 = ln(A+). Another key property is that 
C is a mixing matrix, this means that there exists a decorrelation time a G N such that C" has all entries strictly 
positive. This means that we can connect any two element of the Markov partition in a time steps. 

For every point s = {sjj^g^d £ {l,...,n}^ we can consider the Cartesian product Qs = x^^jdQg^ C T/v. It 

is clear that the collection of Qs with s e {!,..., n}^ forms a Markov partition for Sq. Note that it is natural to 
index the element of this partition with the element of {!,... ,n}^ so that we can associate to this partition the 

symbolic dynamics cq : Z'' x Z = Wf^^ Tpf naturally defined from Cq. We can still call C the compatibility 
matrix and {!,..., n}p the set of possible sequences (namely a = {(j^,i}^^^d is in {!,..., n}p if and only if 

^o-{,io-^,j+i = 1 for every ^ G Z'' and i € Z). Given any point (^,«) G Z'^^^ we will call ^ its space component and i its 
time component. 

The key observation is now that the sets h^{Qs) form a Markov partition for Se- This implies that the space of 
symbolic sequences for Se is the same as that for and that the symbolic dynamics Cg for Se is given by Ce{a) = 
he{co{(T)). Clearly is still Holder continuous. This completes the construction of the Markov partition for Se- We 
thus obtained that the manifold T can be mapped to {1, ... , n}^ ^ where d directions of the lattice Z^^"'^ represent the 
d directions of T = (T^)^ and the last represents the time evolution. Indeed the map Se on the space {1, . . . , n\^^ 
becomes the shift on the time direction, to be called r. 



2.3 SRB measure 



Let now consider the SRB mcasiu-e j-iff^ as defined in section 1. In this case we need to keep N finite because it is 
not easy to give a meaning or construct directly the SRB measure for N — oo. 

Let rn^^ be the measure on {1, ... , n}^^^^ defined as m^°(A) = /i^^(c7^(A)), i.e. is the image of fi^'^ 

via symbolic dynamics c^. The measure can be described efficiently through its restrictions to finite subsets of 

Vn X Z. 

Given A c Fat x Z, m^f ^(ctaIctac) will denote the probability of the event {ct'Icta = cta} conditional to the event 
{(j'\(j'\a = tJAc} w.r.t. to the probability measure m^^, where A'' = (Vjv x Z) \ A and cta is the collection of the a^^i 
for (^, i) G A. 

^From the theory of SRB measures, see [Si] , [GBG03] , it follows that is a Gibbs measure and its conditional 

probabilities satisiy: 
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{ce(T-Krj"))_ 



(3) 



where a' (resp. ct") is the configuration coinciding with cr^ (resp. a'^) on A and with a\c on A'^; t is the image of 
Se through Ce {i.e. it is the one step shift in time direction); Ve^"\ip) measures the expansion of the volume on the 
unstable manifold at the point ip. To be more precise let {wi%iip)}^eVN be a basis on E^^. We will construct one 



such a basis in section 4. Then we have 



def 



\ 



det^r, [{DS^ ■ {DS^ <| 



det|„ [' 



(€) 



(4) 



e,+i 



where u ■ v represent the usual scalar product in M " and det^rj is the determinant of the expression in square brackets 
thought as a matrix indexed by ^ and rj. 

Using the invariance of under Se and introducing the unstable Lyapunov matrix C{ip) satisfying the equation 

DSeWwi%{i^) = 
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we can rewrite the above expression as 

Now the first ratio in eq.(5), when inserted in eq.(4), is vanishing ; indeed the uniform Holder continuity of 'w^^\_{he{4>)) 
and the fact that a' and a" are asymptotically identical in the past and in the future implies that: 

^lim^ (in V'dct,,(u;(«| • wi%)) {c^ir^ a')) - ( In ^dct,, (u;^^ ' ^'i?!)) (ce(r^a")) - ; (6) 

thus the choice of the basis in does not change the result, namely the SRB measure does not depend on the choice 
of the metric as is to be expected from its definition. Calling A^{tp)'''= (^lnC{he{ip))j we finally get: 



A«(co(t^^t')) -A?(co(rV")) 



(7) 



4-00 

E E 

j=-oo few 

Here we used the the fact that = o cq. Furthermore the Holder continuity of A^(co(cr')) implies absolute 
convergence of the sum in eq. (7) because only points asymptotically equal both in the past and in the future are 
compared. 

The crucial point of this construction is that the matrix L(?/))'^=^£(/ij((,y)) is analytic in s due to the fact that 

depends only on w^g^\_{hs{^))- As we already noted w^f'!^{he{'ip)) are analytic in e. We will prove this fact in section 4. 

In section 6 we will apply to eq. (7) the standard methods developed in the study of Gibbs measure in Statistical 
Mechanics. To do this we will need to decompose the "interaction" A^(co(cr)) as the sum of potentials depending only 

on ax"^— {crj}j(^x where X is a finite subset of Z*^^^. More precisely we will decompose: 

(l,j)eVwxZ XcVnxZ 

(these two series are not convergent: they represent the formal expression for the "Hamiltonian" of a Gibbs measure. 
See section (5.2) for a more precise statement). We shall show that we can choose (f)x analytic in e, translationally 
invariant in space and time directions and decaying exponentially in the tree distance of the set X, namely the length 
of the shortest tree connecting all the lattice points in X. In this way (7) can be written as 

gS|^=exp{- E ['Pxia'x)-<Pxia'^)]}, (9) 



SO that one can finally write: 



TOjv (ctaIctac) = — — ^= — . (10) 



This will allow us to show our analyticity claim uniformly in N. 



3 Perturbative Construction of the SRB measure 

In this section we construct the conjugation and prove that it is analytic in e. The technique we use consists in 
expanding he as a power series in e and writing a recursive relation linking the n-th order coefficient to the coefficients 
of order i with i < n. This naturally leads to a tree expansion of the form usual in perturbation theory for Quantum 
Field Theory, i.e. the tree we will introduce are the "Feynmann graphs" of our theory. See also [GBG03] and reference 
therein for similar application to KAM theory. 
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3.1 The Conjugation 

^From now on we will identify functions from Tn to Tjv with their lifts as functions from R^^" to R^^" . Using the 
definition (1) and looking for he of the form /le(^) = ijj + (5/ie(^), we find: 

She oSo-Sqo She = efo (Id + Sh^) . (1) 

where Id is the identity map. 

Setting A = A_ = A^^ and writing /(V') = ^ Z^" ('^^)^ofa '^^^ similarly for 5h^ , we get: 



(2) 



' (S'oV')-A-i(5/if (^) =e/«+(v^ + <5/ie(V')) 

< 

8hi~ {Solp) - \5hf (^) = + ^/le(^)) 

Both equations can be implicitly solved by iteration: 

Shfii,) = -a£^A^'+'"«/«"(5o"(^+^-''"V + '5/ie(5o^^+'"''"V)) , (3) 
p>o 

where Pa = (1 + a)/2. 

It is easy to see that the series in eq.(3) is absolutely convergent, since A < 1 and / is bounded. Expanding 
f^" [ip + 6he{ip)) in power of its argument we find: 

fc>l S=l fcl,...,fc3 



where we have introduced the index x = {£,,a), with a = ±, and J^xi...x,, _ . . •9x,.f with the partial 

derivative in the direction of Wq^^. Moreover we use the convention of summing on twice repeated indexes. The first 
order coefficient of the expansion of the conjugation is then: 



= (-a) E A^+''°r(5o"^^+'-''°V) , (5) 

p>0 



while the fc-th, k > 1, coefficient turns out to be: 

k ki + ...ks=k—l 



shu^)=Y: E (-")E^'^''"(^^-^^'^rfc\)---^'^f4))(^o^'^'"''"^^) • (6) 

S=l ki,...,ks p>0 

^Prom eq.(5) we see that She is in general non differentiable with respect to 'ip. Indeed already differentiating 

with respect to we find a non converging series. On the contrary it is clear that eq.(5) is Holder continuous in ip for 

every exponent /? < 1. 

We can interpret eq.(6) graphically as shown in Fig. 1. 

The l.h.s. of the graphical equation in Figure 1 represents while the r.h.s., representing the sum in eq. 

(6), is a "simple tree" consisting of a "root" r, a "root branch" Xy = {r,v) coming from the "node" (or "vertex") v 
and Sy branches "entering v" , to be called A„^ = (v,Vi), i = 1, . . . ,Sy. 

Even if the drawing in the figure does not carry them explicitly, we imagine that some labels are affixed to the 
node v: more precisely x(u) = {(_{v),a{v)) 6 V/v x {±} and p„ e Z+. Furthermore a label = {^x,ctx) € Vn x {±} 
is associated to each branch A. In the figure above xa„ = x and xa„, = Xj, z = 1, . . . , s. 
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Figure 1: Graphical interpretation of (6). 



The node v symbolizes the tensor with entries 

where p{v) = a{v){py + 1 — Poi(v))- Observe that, in order for eq. (7) and (6) to be non zero, we must have 
ICa„ — CC^')! < 1) due to our definition of the coupling /. 

The line At, exiting vertex v symbolizes the propagator, that is simply (5x;^^.x(t))- 

The line with label x exiting from the bullet of the l.h.s. with label (fc) represents Sh^i.^{ijj); the branches with 
labels Xj exiting from the bullets of the r.h.s. with label (ki) represent 5/i^^^(5'o'"^V)- 

Even if it is not explicitly written in the figure above, a summation over the free indices x{v),x\^. has to be 
performed (note that the summation over x(w) simply fixes x(w) = x, because of the presence of the propagator 

<5x(j;),x)- 

Since eq. (6) is multilinear in Sh^% we can just replace each of the branches exiting from a bullet with the same 
graphical expression in the r.h.s. of the above figure. And so on, until the labels (fc^) on the bullets {top nodes) become 
equal to 1. In this case the end-points represent that is a known expression, see eq. (5), and we will draw these 

known end-points as small dots. 

Thus we have represented our Sh^,.-^ as a "sum over trees" with k branches and k nodes (we shall not regard the 
root as a node) of suitable tree values. In Fig. 2 we draw a typical tree 9 we get via such procedure. 




Figure 2: A tree 9 of order A; = 10 appearing in the expansion for Sh^. Labels ^{vi), a{vi) and p^. are associated to 
all vertices Vi. 
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Note that a tree induces a partial ordering among its nodes: a node w precedes v (and it will be written w < v) if 
there is a path of branches connecting w and v with the arrows pointing from w to v. 

Let us now summarize the discussion above. Let 7fc(x) be the set of rooted trees with k branches and k nodes, 
with labels x{v),pv attached to their vertices and x{vo) = x, where vo is the last vertex preceding the root. Given 
9 e 'Jfc(x), let the value of 9 be defined as 



„/ \ fx{v), x{vi)... x{vs^) 

Val(^, V) = n ( - "(^)) n (^o' V) ' (8) 



fx{v), x(di)... x{vs^) 

where V\,...,Vs^ are the nodes immediately preceding v and p{v) = J2w>v + 1 ~ Pa{v))- With these definitions 

Sh^f,-^{tp) can be calculated as Sh^j^-^^ij)) = X^^^^j^^^^ Val(^, ■;/')• 

3.2 Convergence and Regularity of the Perturbative Expansion of the Conjugation 

By definition depends only on tpnn{0) so that it is analytic in V'^={tp^^ e C| |Im^^| < Tq, i = 1,2, ^ G nn(0)} for 
some ro > 0. Calling G the maximum of g on T>, from Cauchy's formula we get 

\r -^--^m < G^^^i-^ < g4. (9) 

^0 ^0 

where mi, ...niD are the multiplicities of the partial derivatives with respect to the D'^=2{2d + 1) = 2|nn(0)| possible 
variables, (thus mi -{-... + mo ~ s). 

In the same way, if ij) and V'' are identical on each site but ^' € nn(^) and if < /? < 1 we get 

ir--- -^(5?V) - r'^-^^(5oV)l < G^^±ll:(27r2)^A-''^|V'e - (10) 

^0 

where we have used the periodicity of /. Next we bound the value of a tree 9 G ^^(x). Using eq. (9), for 6 G '?fc(x), 
we find: 

||W(0, OIloo < n 5r = -£r n ' (11) 

where we used that, if ^ G '7fe(x), X^„gg s„ = — 1. 

The sum over the trees can be interpreted as a sum over the topological trees and a sum over the labels attached 
to the trees. If 6^ is the set of topological trees of order k, we get: 



ii^/^(.)iu<E E ^En^-^E E ;^(t^)'=^ (12) 

< E 2'=(2d+l)'^^(^)'=<22'=2^(2d+l)'=-^(^ 



where we used that: 

1. 2*^ is the number of terms in the sum over the a{v) indices; 

2. {2d + 1)*^ is a bound on the number of terms in the sum over the values of the ^(v) indices not making Val(0, ip) 
vanish (observe that, given a tree 9, its value Val(^, V') is vanishing unless \$,{v') —5,{v)\ < 1, where v' is the node 
immediately preceding v); 

3. 2^*^ is a bound on the number of unlabeled rooted trees with k nodes. 
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In the same way we find that, if and are identical on each site but ^' and if < /? < 1: 

^ E (t3^)' (2^')'^^2'=(2d+ + 1) < (13) 

^ (t^It^^)' (2^')'^^2'=(2d+ l)'=(2fc - 1) . 

So the map /le : Tjv — * 7iv exists, it is holder continuous w.r.t. V and analytic w.r.t. e in the complex disc |£| < 
with ^ 

(14) 

In order to prove that he{tp) is an homeomorphism. wc have to show that it is invertible. The proof is easy and 
standard. Regarding injectivity, note that two distinct points i/ii, -02, are necessarily far order one in the "future" or in 

the "past", namely there exists an integer n S Z such that I-S'q^/'i — SQip2\ = 0{1); then S" (jie{'4'i)j ~ (j>-s{'4'2)j = 
Sq{iPi — -02) + 5he{SQtpi) — 5h^{SQip2)^ cannot vanish as the first term is order one, the other two of order e; thus 
it cannot be but he{ipi) ^ he{^2)- Regarding surjectivity, since / is a continuous injective mapping on a torus, / is 
necessarily surjective (the proof is trivial on T and it can be easily extended by induction to Tjv). 



£/3 = 



4 The Unstable Direction 

In order to explicitly compute the SRB measure, we have to construct a basis for the unstable subspace , and the 

expansion coefRcient V^'' associated to it, as explained in section 2.2 above. Note that we cannot use to find a 
basis for because it is only Holder continuous. 

To find the unstable base {w^^]_(V')}{ey and the Lyapunov matrix C{ip) we have to solve the following equation 

{ds, (V) = 4fi {Se (V')) -cf' (V) . (1) 

In general this equation cannot have solutions analytic in £. In fact, from the general theory we know that the unstable 

vectors {7«^^^('0)}^gy are not differentiable with respect to 0. But, as we previously pointed out, to compute the 
SRB measure we need only to know the expansion coefficient at the point he{ip), i.e. 'Di"\hs{'4'))- Let us define 
wl^l_{he{ip))'^= vi^\ip) for ^ e y and note that vi^\ip) satisfies the equation 

{DS,){K{i;))vi^\ilj)=vii^{Soi^)Li^ii^) , L(V) = £(/i.(V)) • (2) 

We will show that this equation admit a solution analytic in s. Moreover the determinant of L{tp) is all what we need 

to compute the SRB measure. 

At this point, it is convenient to write eq.(2) in components. Denoting by y the double index rj^ (again x'='^^"), 
defining 4^^(V')'= K*,i'(V')4?/3 and {DS, 4'/^) (V--) =^ Ex get 

sf^y (/le(V)) K^,^; W = v}^^ {Som^'W ■ (3) 

Now, defining the corrections 6L and SV as follows: 

L^\^pfiJ'x-^6^, + SLi\i^) , K^^J{^pf^V,^^2+SV^i\^P) with V^^^^=6^,r, , V^f^. = , (4) 
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we find that (3) is equivalent to: 



+ sf^^'y{h,{i,))6Vj'\^lj) - 6V^i\Soi')SLip{i,) 



(5) 



Of course tiie above equations cannot determine completely the basis and its associated matrix: indeed given a solution 
{Vy\i;)}, {L^Pi^)} of eq.(2) and a generic invertible Holder continuous matrbc R^p{iI>), also {vj"'\il>)R^P{il>)}, 
{R-^'^^{SQ'il>)L^^{^p)WP{i))} solve (2). Thus it is possible to add a constrain to SVj^\^py. a possible choice, which 
simplifies a lot the expressions above, consists in taking 6V''l\ip) = 0, so that (5) becomes: 



+ 6V^^9iSoij)6Lip{ij) . 

An implicit solution of (6) (to be inverted iteratively by a new tree expansion, see below) is: 
' 6LiP{iP) =efi^^P^{h,m+ef^^'^-{h,m6V^'^\i,) 



(6) 



(7) 



As for the construction of the conjugation, we can expand in power series of e both sides of eq.(7) and equate the 
coefficients of the same order, thus finding an iterative solution of 5L(fe) and (5V(fc) . The first order coefficients are 
given by: 



5L\Pp)=f^*^p\,l^) 

j>0 



while, for A; + 1 > 2: 



s>l, ki>l 



(8) 



(9) 



and 



+ 



+ 



fciH i-ks=k 



feiH \-ks=k 

E ( 

s>l, ki>l 



,r^SV^9,, .dh"^? Sh'^C ^) (V') 

ri (fei) (fes) (k,))\l'' 



E Ea^^'"^( 

L, ki>l i>0 

■■+k s=k fS^.v 

E E^^^"^(^ 



+ 



s>l, ki>l j>0 



fg — lV V (kl) (fe2) {ks)J^ U V/ 



+ 



(10) 
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(s - 1)! 

s>l, fej>l j>0 ^ ' 

s>2, fej>l j>0 ^ '' 

This two relations, together with (6), allow a recursive construction of and SV . Obviously, repeating the discussion 

of §3.1, one finds that bL and SV can be expressed as sums over trees, obtained by suitably modifying the construction 
of previous section. It can be easily realized that the estimates for the tree values are qualitatively the same as before, 
see eq.(ll),(12),(13). We point out the differences appearing in the tree expansion for SV: 

1. the nodes can be of 4 different types (corresponding to the 4 lines in eq. (10)); so that the number of possible 
labels for a tree of order k is larger of a factor 4*^; 

2. the number D„ of derivatives acting on a node function can be either s„ or s„ + 1, see eq.(9),(10), so that 
differs from the combinatorial factor St,! by at most + 1; then the final estimate contains a factor that can be 
bounded by (1/ro) Ht-Csf + 1) < e'^/^'o- 

A similar discussion can be made for the tree expansion of bL. 

The result is that L and V are analytic in e and Holder continuous in with exponent < /? < 1 in a disc |£| < 

with e^Jj smaller than the convergence radius e/3 of h^, see eq.(14). Note that also in this case e^Jj is independent of TV. 
As already explained (see section 2.2 and in particular eq.(7)), in order to compute the SRB measure we need 

A« = (logL)« =^ - log A + where: 

(5A«(V')= [log(J + AJL)1^^ = y ybL^'^H^^-'-bU^'-^^iM,) . (11) 

L J ^-^ s 

S>1 

(no summation on ^ is intended). Expanding eq.(ll) in series of £, we get: 

S^UW = ''^^I1 (12) 

s>l, ki>l 

Again, the last equation, together with (10) and (6), allows a recursive construction of the coefficients SA^^j^-^ and the 
result is that A^ is a sum over (suitably modified) trees. The bounds are still qualitatively the same, so that A^ is 
analytic w.r.t. £ in a suitably small complex disc (independent of N) and Holder continuous w.r.t. ip. 



5 SRB Potentials 

The next step towards the constriiction of the SRB measure and the proof of its analyticity consists in the expansion 
of A^ in potentials (px- From the analysis of previous sections follows that A^, as well as h^, V^^^ and L^^, can be 
expanded in convergent sums over tree values. We will discuss here how to expand h in potentials, since the analogous 
expansion for V, L and A is conceptually similar, just more involved due to the more complex structure of the trees. 

We will proceed as follows. We first write the values of the trees in terms of the symbolic variables a. We then 
decompose each of these values as a sum of terms only depending on the a's on finite but arbitrary large sets. Finally 
we define the associated potentials by collecting together the contributions which depend on the same cr's. Our goal 
is to obtain potentials defined over sets with rather arbitrary shape but decaying exponentially with the tree distance 
(see after eq.(8) for a precise definition) of their support. 

To begin with we expand the derivatives of the perturbation function / via a telescopic sum. Given the digits s 
and s' € {1, . . . ,n} we can always find a sequence of digits S(s, s') = S1S2 ■ ■ ■ Sa-i such that the sequence sE(s, s')s' 
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compatible, i.e. such that Cs^^si+i = 1 for z = 0, . . . o — 1, where sq = s and Sa = s' . Chosen a sequence a € [1, ... , n}^ 
once and for all, given a G {1, . . . we can define its restriction to time j, as follows: o"| ^ = a^ i if \t\ < j, 

o'lt = CT^.t if \t\ > j + a and the gap is filled with the sequence constructed above for s = cr^,±j and s' = <5'5,±(j+a)- 
We can now define: 

=^ X] (^"""^ (€) ) (1) 



where ^ is the spatial coordinate associated to X and nn^-'^(^) =nn{^)xlj,lj = [— j,j]nZ. Since |co(cr-') — co(<T-' "'^)| < 
'0) 



cA-' for some c > 0, f^-^^ ^' is bounded by: 



||/(7)^^-^1U<G^^cA^-. (2) 
5.1 Decay of the Potentials for the Conjugation 

Inserting expansion (1) in the definition of the value of a tree eq. (8), we find: 

px(-u), x(i)i)... x(i;s^) 

Vai(^,co(a)) = n E (-a(^;))A^-+^°('-) ^^"^ ^ , (t^(''V„„o„,(j(,))) , (3) 

t;eej"„>o ^' 

where we recah that p„(„) = {l + a{v))/2, p{v) =J2w>v°^i'")(.Pv + ^-Pcx{v)) and nn'^^\£,) =nn{^)xlj,lj = [-j,j]r\Z. 
The above expression can be seen as a sum over the values of a new kind of trees, identical to the ones described in 
section. 3.1, but with a new label jv e N attached to each node. Let 7fc(x) be the set of these new trees of order k 
contributing to Sh^j^^ , i.e. 6 G Tk (x) is a tree with k branches and k nodes (the root is not a node) with the following 

labels attached to the nodes v € 9: ^{v) €V,pyG N, jy € N and a{v) G {—1, +1}- 
Given 6 GTk{x), its value is given by: 

„x(t;), x(ui)... x(t)s„) 

Val e{a) = n ( - c^W)A^"+^°"" ^'"^ ^ , (t^^''V„„(,„,(5(„))) . (4) 

SO that we have 

<5/if,)(co(a))= ^ Val0(a) with ||Val 0||oo < (^) ' H • (5) 

where we have used eq.(2). We can now define the support X{9) C V x Z of a tree 9 £ ^(x), as the support of the 
spin variables on which Val0 depends in a non trivial way, plus a center 0). More precisely: 

X{9P{{^,0)}u[jc{^{v),p{v),j,) , where c{^{v),p{v),j,y^^ [j [j {v,p{v) + i) . (6) 

veO rienn{${v)) \i\<jv 

namely C{£^,p,i) is a cylinder centered in (C,]?), with the spatial base equal to the set of nearest neighbors of ^ and 
with height equal to 2j. Then X(9) is the union of (^, 0) and of cylinders of this kind, one for each node v of the tree. 
The point 0) has the role of center of X{9) and is added to X{9) for later convenience (note in fact that Val 6{(j) 
could not depend on cr(j,o))- 

Given a set X C if^^ we can partition it in a natural way as a union of time-like segments. More precisely given 
^ e Z let = {(5, i) G Z^ I i e Z}. The intersection between and X can be uniquely partitioned as a union of 
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maximal connected segments. The collection of all this segments forms a partition of X in nx time-like segments 
{Ri{X)}i^i nx- Let now the center of Ri{X). If F is a subset of if''^^ we call tree distance of Y, dt{Y), the 
length of the minimal tree connection all the points of Y . Finally let dc{X) the tree distance of the set {ri]i=i^,,,^nx- 

^From the previous bound on the value of a tree 6 G 7fe(x) eq. (5) can be interpreted as the tree distance decay of 
the contribution of order k to 5h. Indeed: 



kl'l|ValO|U<(;2£N!f)' 



{xm\ 



(7) 



where: 



1. the factor A'''^*^''^*^^'' comes from Hree ^^^''"S in f^ct py is the displacement in the time direction of the cylinder 
associated to the node v w.r.t. the one associated to the node v' immediately following v, and 1 is their maximum 
displacement in spatial direction, so that X]^(l +Pv) > dc{X{6)); 

2. we used that nx < (2rf+ l)k in order to bound \e\^/'^ with lel^c^d+i) . 

3. the factor aI^'(^W)I comes from XI^^qX^-; 

4. the global power l/(2d+ 1) in (7) comes from the size of the base of each cylinder, namely we used the fact that 
the number Ux of segments is less than 2d+l times the number of cylinders in X{6). 

Collecting together all the trees 6 which have support X{6) = X for a given X, we get: 

X(0)=X 

5h^{co{a))= J2 ^f'xM with 6h^{axf=J2^'' E Val ^((t) . (8) 

X9({,0) fe>l e<£Tk{x) 

So, using the bound (7) for |£| small enough, 70 = 2(23+17' '^0 ~ ~27ologA, vq = |e|'''° and a suitable c > we get 

nx 

||5/i$:||oo < ce-««'^=(^V^'^ ]Je-'=°l-^*WI , (9) 

i=l 

namely 6hx decays exponentially with the tree distance of X. 
5.2 SRB Potentials and Their Decay 

Proceeding as above for the function A^(co(cr)) we obtain that we can write it as: 

Xc{Vnx1') 

where by construction (/'^''^^ is different from only if (^, 0) G X. The function (f>x'^^ is again given by a tree expansion 
analogous to that in eq. (8). Moreover we will set 

<Px i'^x) = <|)'^-/x{(^x)■ 
We can define: 

{U)ex 



14 



so that we formally obtain eq.(8), namely, given It = [—T/2, T/2] n Z (T even) and calling A = VJv x It, 

^ A«(co(tV))- J2 M<^x)^0{dA), 

where dA is the boundary of A and the correction can be exactly computed from the definitions above. 

Note the potential 4'x{o'x) is invariant under time and space translations (respectively for the definition of 
and for the periodic boundary conditions), namely: 

(l>x{(Tx) = (l>pSrjx{crx) for any j) e Vjv x Z. (10) 

Moreover it can be bounded by: 

nx 

||<^x||cx> <ce-'^^'^=Wz/I'^ J]e-'^il«^WI , (11) 

i=l 

for suitable c, 71, ki > and vi = \s\'^'^. 

6 Analyticity of SRB measure 

In the previous sections, we wrote the SRB measure as a Gibbs measure with translationally invariant potentials (j)x , 
decaying as in (11), and with hard core interaction in time direction. Moreover the potential ^x is analytic in e in 
a small disc in C around the origin (independent of A''). A well known technique to show analyticity of the Gibbs 
measure w.r.t. e is the so called cluster expansion. 

U A=Vn X It, with It'^=[-T/2,T/2] n Z for some even T e N, we call Aa = Vn x lT+2a- Given a boundary 
condition &£{!,..., , we define the pressure Pa as 

PA=Vr'l0g^e"^-nA^«'^^^"^^ . (1) 
cr 

where the sum is over all the a that coincide with cta on A, to a on A^ and with T,{a^^T/2,^^,T/2+a) in the space 
remaining. It is well known that the pressure Pa can be considered as the generating functional for the Gibbs states. 
From its analyticity our main theorem will follow easily, as we will see in section 6.4. 

6.1 Decimation 

In presence of hard cores we cannot proceed in the standard way (Mayer's expansion), since the standard proof, see 

[GMM73], requires weakness of the original interactions. We can overcome this obstacle by a decimation, see [C081], 
namely considering the statistical system on scales larger than the length of decorrelation of the hard core. 

6.1.1 Decimated Lattice Ad 

For each ^ e Vjv, we divide the time interval It'^={Q x It into an alternating sequence of blocks, called of "S-type" 
and "if -type", B^\h^^\ B^^\h^^\ .. ., B^^~^\ H^~^\ B^^\ containing a number of spins respectively equal to 6 = 1 
and h = h^a — 1, with /iq G N to be chosen later. For this reason we choose the number of points in 7|, to be 
= thoa + 1, namely T = ihoa. 

Remark The choice = 1 is special for the present case, in which the unperturbed potential is vanishing. In general 
one could treat with the same technique the case in which the unperturbed potential is order one, with a sufficiently 



15 



fast decay of the tails, and in that case b should be chosen suitably large, see [C081]. Such a case arises, for instance, 
when the unperturbed system is the product of non linear Anosov maps on T , namely in the case treated in Appendix 
A. The present discussion could be easily adapted to cover that case. 

Let /3f ="'^o-(| _T/2+i/ioa), ? e Vjv, i = 0, be the spin in the block fif and =^{<^(€,p)}(j,p)eHW, ^ e Vn, 

z = 0, 1, be the collection of spins belonging to the block ijf ; it will be regarded as a sequence of h P spins: 

= (Pi{ri^^), . .. , I3h{r]^^^))- The lattice obtained considering the H and B blocks as points: 

will be called the decimated lattice] on A^ the distances will be computed by thinking of it as having its sites spaced 
by 1 also in the time direction. 

If X C A, Y{X) will denote the corresponding subset in Ajj, namely the smaller subset Y C A/j such that the 

union of the B and ff-blocks in Y contains the set X. Defining $v'(/3y-, Tyy) = Y^x:Y(x)=y ^x{ox)-, eq.(l) can be 
rewritten as 

where = C'/3/3i(»7)C/3i (,,)/32 (»))••• C'/3h_i(»7)/3ft(,,)C/3,.(,,)/3' • (3) 

Observe that, from eq. (11), if Y do not coincide with a single if-block, $y satisfies a qualitatively equivalent bound: 

UY 

ll^rlloc <ce-'^''»(^)l7"^ JJe-'^l^*^^)! , Y ^ Hf , (4) 

i=l 

for some c, k,7 > and u = \s\^. Whereas if y = ilf for some ^ gV and some i = 0, . . . ,£—l,we have ||$r||(x) < hu. 
6.1.2 Averaging over many degree of freedom: the Perron-Probenius's theorem 

Decimation is a renormalization group technique, consisting in summing first on the if-type spins, thus getting an 
effective statistical system for the B-blocks: the idea is that if the B-blocks are sufficiently far apart, after the averaging 

of the ry's, the /3's should be almost independent, as if there were only small interactions among them. The technical 
tool we shall use to prove rigorously that the effective interactions between the /3's are small is Perron-Frobenius 
theorem. 

Let Z{(3, (3') be defined, with a little abuse of notation, as: 

^(/?,/?')'= E^(/^'^'/5') = C;^? (5) 
n 

Observe that 1 < C^^, < q". Since C has strictly positive entries, we can apply the Perron-Frobenius theorem and 
obtain that C" and its transpose C°"^ admit a non degenerate eigenvalue I > with eigenvectors tt and tt* respectively 
such that TTcr, TT* > for any a = 1,... g, and ^^T^^T^a — 1- The eigenvalue I is maximal in the spectrum of C"; 
namely, if we define P as the projection matrix P^ct' = ^aa' — tTctTt*/ , we have: 

||(Z-iC"')'=Pu;||oo <c„e-"'=||u;||oo , (6) 

for any uj G M*^ and with 

- log(l - [min(C^,„/0)]') > 9"'" • (7) 
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As a consequence, 

Z(/3,/3') = Cfpl = ^C|^(7r.7r^, +P<.^0 = 

(7 

= / °7r/37r^, 1 + — — =/ "7r/37r^,e ^''"^ ^ (8) 

with J(/3,/3') = 0(e-'''"2" °). It is now clear that taking ho big enough we can make the two body potential (3') 
as small as needed. 

Using eq. (8), introducing a new effective potential W including the contributions from $ and 7, defining 

i=0 j=0 * 

and using limA^oo logrif.j S/3('> ''''^^ ■* = (as it follows from the normalization condition T^a^a = 1)) 
we can rewrite Pa as 

PA = ^log/ + i-log^^m(/3A,7?A) n e-^-^^-'"-) 

' ' Pa Va ycA 

with mi(3^,VApl[I[ ^ .-umM li 'JJ.)' (10) 



where m{j3\,r]\) is a probability density. Observe that, if one chooses ho ~ — logi^ (so that both hu and I{f3,f3') are 
small), the new interaction W satisfies a bound similar the one of 



||Wy||oo <ce-«'*=(^)F"^[]e-^l«*(^)l , VFcAb, (11) 

for some c,k, 7 > 0, F = jel'''. 

6.2 Mayer's Expansion and Polymer Lattice Gas 

We shall now expand the small potential appearing in the expression for Pa, via a Mayer's expansion, obtaining the 
pressure for £ = plus a correction. 

It will be convenient to collect together the contributions of the potentials whose supports have the same closure, in 

the following sense: for a set formed by a unique point H^^^ G An we define its closure as (//f ) =^ (Pf , ijf , 5**+^^) 
while for a set formed by a unique point P^'^ G An we define {B^^)'''= (P^*^); finally for Y c An we define its closure 
as Ft/ UGevJC). 



We say that a collection C = {Ym}m=i of sets Yi C 'Zf^'^^ (think to them as molecules) is connected if, given a couple 
(^n, yfin) e C X C, it is possible to find {Y^^ j^^^, such that Fin n Yrm ^ 0, F,„, n F,„,^, ^ and F™^ n Ffin 7^ 0. 

Writing e"^'''^^^'''^-' as the value for £ = plus the correction, namely 1 + (^e~^^ ^^^^ '"^^^ — Ij, expanding the 
product over Y c and collecting together the connected components, we can rewrite eq.(lO) as: 

J'A-^log/ = ^log^^m(/3A,r;A) ^{^)\{p{l\P^,ri^) , (12) 
' ' 13a VA rcAo 7er 

where: 
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1. 7 is a subset of Ajj, to be called in the following polymer (they are, indeed, the union of a connected collection 

of molecules); 

2. r is a collection of polymers: F = (71, ... , 7„), n > 1 and F C means that 7 C Ad, V7 € F; 

3. T(r) is the function equal to 1 if 7 fl 7' = for every 7, 7' e X with 7^7' and otherwise; 

4. p(7|/3^, r/.y), 7 c A£) is defined as 



P{1\P.,V,)=' Eli E nl-"^^-^'^""^^^-!) (13) 
q>i ^' n^-.y? i=l 

where the * on the sum means that Fi, . . . , 1^ is a connected collection of subsets of A^; 
5. the term corresponding to F = must be interpreted as equal to 1. 

The key observation is that, thanks to the above definition of closure, in (12) we can sum over rj spins before 
summing over the P spins. After doing this the measure m(/3A,??A) factorizes, i.e. 



def 1 



def 



^ ^ m{Py, ?7^)p(7|/?7' 



Namely we have rewritten Pa as the pressure for £ = plus a correction having the form of the pressure of a "polymer 
lattice gas", with activities ^(7) and hard core potentials T(F). 

6.3 Cluster Expansion and Its Convergence 

A standard argument, exposed for instance in [GMM73], [PS99] or [GBG03], leads to 

Pa - i log / = log ^ T(r)p(F) = ^ ^ ^^i^)Pi^) ' (15) 

where is the Mayer function, defined as 

^ {li,---,ln)= { , (16) 

[1 if n = l. 

where G{n) is the set of connected graphs which can be drawn on n vertices labeled 1, . . . , n by connecting with links 
couples of distinct vertices; the function fi'Jij'yj) is equal to 1 if 7^ H 7j 7^ and otherwise. By construction, T^(F) 
is different from zero only if F is a connected collection of polymers. Observe that F could contain many copies of 
the same 7. More precisely here F represents a function from the subsets of A£) to N (and we can think F(7) as 
representing the number of copies of 7) such that ^("f) — 

A bound for ^(7) can be obtained as follows: 

K7)l<l|p(7h-)lloo<E^ E nir^-Jlooell^-^ll-, (17) 

P>1 Yi, Ui'Kj=7 *=1 
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Using the bound (11) (and that, if U^^^F, = 7, one has J2i=i W^yAI < ci'\'y\) we find: 



P- — 



(18) 



i=i 



We can now use the connectedness constraint on the sum in order to extract a factor exponentially small in the size 
of 7. Indeed, if U^^^F, = 7, one has both d'c{Yi) > dc{Yi) and ''^y, > n-j- After extracting such a factor we 

can relax the constraints on the sum, so that: 



|p(7)|<ce-l'^l 



\\R^m 



It is easy to see that the last sum is bounded by c\-)\U^^^, so that: 



g-^|fl(Y)iy 



(19) 



|p(7)| < ce^'^l'^l 



Nc(7),-7n^/2 



i=l 



(20) 



p>i ■ 



i=l 



for some c, K',7' > and = \e\'^ . Using the preceding bound we can easily prove that: 



sup Y. |P(7)I <c(i^')'^'e-('^'/2)r _ 



(21) 



where ^(7) is the diameter of the polymer 7. A standard theorem, proved for instance in [GMM73] [GBG03] , states 
that, if p{"f) satisfies (21), then 

sup T^^(r)|p(r)| < c{v'f'^e-^^''^'^'' . (22) 

5{T)>r 

This implies that, varying A, Pa is a uniformly convergent sequence of analytic functions in a domain independent 
from A. The limit, still analytic in the same domain (thanks to Vitali's convergence theorem), is independent of 
the way the thermodynamic limit is performed {i.e. one can send the time side of A to 00 either before the spatial 
side is sent to 00 or together with it), thanks to the exponentially fast convergence of the sequence, implied by (22). 
For the same reason, the limit is also independent of the choice of boundary conditions and, because of translational 
invariance, it is equal to: 

r9(o,o) 



Pt/ hm P^ = -\ogl+-^ 

|A|^cx) a haa 



E 



(r)p(r) 



(23) 



def 

where |r| = | U^er 7I and 2/(/ioa) = lim|A|^cx> |Ar.|/|A|. 



6.4 Analyticity of the Mean Values 

The analyticity for the mean value of an analytic local observable C'(V'y) (depending on the variables in the finite set 
y cZ") is an easy corollary of the previous result. 

We first observe that ^^^^{O) = HmAr.T— 00 i)eiVNXiT) o o SI). This is true thanks to the 

time and space translation invariance of p^^^ . Moreover it is possible to decompose O as 

0(/le(co(a))y) = Y <^x'H^x) 

xn(yx{o})=0 



19 



This can be done expanding 0{h^) in power of ip, using the representation of given in section 3 and 5 and collecting 
the terms with the same support. Moreover we will set: 

and 

pivxnynx^n 

It is easy to realize that Ox is invariant under space and time translations, and satisfies 

nx 

llOxlloo <cyi/"^e-'"^=W]Je-''l-^*WI , (24) 

i=l 

for some k, 7 > 0, u = \e\'^ and some constant cy > which depends on the size of V. Setting A = Vn x It, the 
thermodynamic limit of the mean value of O(V'y) can be written as: 



1 y- e-S^nA^«'^^^"^'-^°^^"^^l 

^sRB(0^ = lim —a^log^'^^ 



Via a new cluster expansion we find: 



=^dMO ■ (25) 



M^--(O) = hm i^ac Yl T^(r)(p^(r)-p(r)) (26) 



where p^ij) are the activities corresponding to the potential (j)x — (Ox- For |C| small enough, the potential (px — C^x 
satisfies the same bounds of 0x so that J2rnAmt$ ^^{^){p''{^)~P{^)) ^ uniformly convergent sequence of functions, 
analytic in e and C in the product of two small discs. This implies that /i^^^(O) is analytic in s and given by: 



r3(o,o) ^ 
p^^d+i I I 
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A Unperturbed Nonlinear Dynamics 

The result about analyticity can be extended to the case in which the unperturbed dynamic is made up of independent 
nonlinear analytic Anosov systems sq : T — > T . We suppose that there exist u±(V') and X±{ip) such that: 

{DsoV+){cP) = X+i(l>)v+{so{4>)) , {D,,v_){cP)=X_{cf>)v^{so{4>)) , (1) 

with (j) s T^, v±{(j)) and A±(0) Holder continuous and |A+((/))|~^, |A_(0)| < A < 1. Then we consider again a 
perturbation /(V') on Tjv analytic in ip. Observe however that in this case the most naive example of perturbation, 

with f~^{il>) analytic, is no longer an analytic perturbation. 
A.l Conjugation 

The constitutive equation for She, hfted on M^^^ , is: 

So{he{^)) +e/(/ie(V)) = So{i^) + Sh4So{^)). (2) 
In order to exploit the hyperbolicity, it is convenient to arrange the terms as follows: 

{DSo 6he){iJ)-5h4So{^|J)) = -ef{iP + She{i^)) - [So{he{ij)) ~ So{i^) - (DSq ^MW]- (3) 

Define 

C(,l,... (71,8 V / Cl = --- Cs=U 



whereas 

X X 

Writing <5/ie(^)'= Ex^/^^C^)^^^, and {DSo «'^S)(V')'= Ey ^o'"(^)"'S(^o(^)), with S^'^'i^) = XU^i)S^,y, we 
get: 

X^Sh-iij) - 5h-{SoW) =-eJ2 C^'Z ^^ ^^e' ■ ■ ■ Sh-^) W - E { ^\r' sh-^ - ^K') (i;), (4) 

s>0 ■ s>2 

Finally, the recursive equation for the Taylor coefficients of Shf^ (ip) is: 

„>0 m=0 s>0 1=1+ ■■+'=8=1= 

+E( n KHsom)^: e (^v-^'^r^o---^'^r4))(^ow^ • 

p>0 m=0 s>2 lci + - ■+fea=fe+i 

l'i>i 

A similar equation holds for x = ^~ . 

^From now on, the construction of the conjugation function goes on as in the Uncar case with similar considerations. 
We have only to take in account the fact that a tree of order k (w.r.t. e) does not necessarily have k branches, because 
of the term on the last line of (5) (to be called a vertex of type 0). Since the number of lines entering a vertex of type 
is > 2, one can easily prove that the number bk of branches of a tree of order k is boimdcd hy k < bk < 2A; — 1, so 
that nothing qualitatively changes in the bounds and the proof of analyticity of She proceed as in sections sections 3 
and 6. 
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A. 2 Unstable Direction 



The perturbed unstable direction in the point hs{ip) is given by the equation: 

Setting 'w^P^{h^{il)))'^= v'i'\il)), it is convenient to rearrange the terms of the equation in the following way: 



Defining vi^^iP) = Ex ^-^,x and using again the considerations of section 4, we finally get 

, yx,yxi---Xa 

s>0 



s>0 

,C'Xy,xi...Xs 



s>l 



(with [DS^ ^°wi^^^{tjj)'^= Sq^'^^'" {tp)wl^'^\So'ip))- Again, because of the third term in the r.h.s. of eq.(7) , the 
number of branches of a tree appearing in the construction of 6V and SL is greater (in general) than the order of the 
tree itself. This is not a problem, since one can easily realize that, again, the number 6^ of branches of a tree of order 
k is such that k < bk < 2k. 



A. 3 SRB Interactions 

Following the proof in section 4 and, proceeding as in sections 5 and 6 one proves analyticity of SRB distribution. In 
fact the only (slight) difference in the construction of SRB potentials is in the telescopic cutting necessary to represent 
h, L, V and A as sums of local functions of spin variables. Notice that now each tree node is associated to the product 
of a node function fv{ip) {e.g. in the case of a tree contributing to 6h, fv can be a derivative of / or a derivative of So, 

see (5)) times a product of local Lyapimov exponents, like the factor 0^=0 -^^^ (''o"(V'c))'=^0+(?'j •^o''"'' ("0?)) in (5); 
the analogous expression appearing in a vertex with a(v) = — will be denoted by Q-{p, SQ^'^^tp^). So the total node 
function associated to a vertex v will be now of the form 

F"(p(t.),5„^('')(^))'=^e„(„)(p(«),.?;(")(V'cJ)/„(5o^("'(^)), (8) 

where u' is the vertex immediately following v. The telescopic expansion (1) has to be done separately for each of the 
factor in the above equation (A('0) is holder continuous), getting in the end potentials with the same kind of decay 
rate. The bounds are not qualitatively changed and the subsequent analysis of section 5 follows so that, by suitably 
modifying the decimation procedure, analyticity of SRB measure can be proved. We point out that a main difference 
in the proof of convergence of the cluster expansion is that now the unperturbed potentials are not vanishing, but have 
support only on time like segments / C Z, and are exponentially decaying with the diameter of /. For this reason one 
cannot proceed exactly as in section 6. The standard way to treat this problem, see [C081], is to choose a length r 
such that the unperturbed interactions on sets I, diam(/) > r, are small enough for the cluster expansion. Then one 
fix the size of the B-blocks b = r, and the size of the ff-blocks, h, such that the Perron-Frobenius theorem is true for 
the reduced partition function Zr{f3^^\ri^\ P^^'^^'^), in which only the interaction on sets / C {B^^^ U H^^^ U B'^^^^^), 
diam(/) < r are taken in account. 
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B Green-Kubo Formula and Large Deviation 



In this section we deal with an appHcation. We introduce the local phase space contraction rate on a volume Vq cVn 
averaged on a time Tq, given by 



,def 1 

I 



|Ao| ^ 



with Ao = VoX lo and Io'^=^[-To/2,To/2] n Z. We prove a Green-Kubo formula for t/Ao, from which it will come 
out that generically its mean value 77+ is strictly negative. Furthermore, we can show the large fluctuations of r]\g 
around r}+ satisfy a large deviation principle, namely they are asymptotically described by a strictly convex free energy 
functional F{r]): it can be obtained as the Legendre transform of the generating functional P{() = P^j^ see eq. 
(25). 

For the rest of the appendix the SRB interaction will be called {^x}xcZ''+^' *° remind that they are derived from 
the unstable restriction of DSe- 

Theorem Bl: Given such that r]+ < 0, 

1. P{C) is analytic and strictly convex in (, for \e\ < eq, \(\ < 1, with £0 small enough; 

2. the Green-Kubo formula is valid: 

a,^P'(0)|e=o = -ia,2p"(0)|,=o. (2) 

Theorem B2: Given such that 77+ < 0, 

1. the free energy F{r]) is analytic in r], for \s\ < eq, and ri € [P'(— 1), P'(l)]; 

2. if[a,b] C [P'(-1),P'(1)], then 

U^log/x"^^(r?Ao e k^>]) = max -AP(r?,r/+) , (3) 

lAoHoo |Ao| ri&[a,b] 

with AF{r],r]+f=^F{r]) - F{r]+). 

B.l Local Phase Space Contraction Rate 

Repeating the construction of SRB potentials leading to (8), we set: 

xnAo^o 

|Ao| 



for a suitable potential (j)x , satisfying 



1 10x1 loo < ce-'^''=(^V"^ JJe-'^l-^^WI , (5) 

i=l 

for some c, k, 7 > and v = \s\^. ^From the invariance under time translations of the SRB measure, we have 

Um 7-— ■ 

IVbHoo |Vbr ' " ^ ' lAoHoo' 



n+ Jim -^/.^--(log det(£>5.K) = , lim m^^^K) 
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lim fj,^^^{(j)x)= (6) 



= lim lim — — - log 



1 V p"SxnA^0'^i('^^)+^ExnAo^0 



lAoHoo lAHoo |Ao| V- p-Z]xnA7S0'^i('^^) 



C=0 



It is easy to show the last expression is equal to the one with the summations over X n A 7^ and X DAq ^0 replaced 
by X C Aq and without the limit in A (since the correction is only a border effect; or simply using again the cluster 
expansion developed in section 6.4). In this way, defining the generating function P{C), as 



P(C) =^ lim ^ log -— . (7) 

^"-^ |AH<^|A| ^ e-SxcA^i(--) ^' 

we finally get: 

77+ = P'(0) . (8) 

Analyticity is achieved by cluster expansion (we do not need ^ small, but we can take, say, \C\ < 1, since {4>x}x are 
0{e)). 

B.2 Green-Kubo Formula 

Consider the case in which so is the Arnold's cat map defined by (1). 

Using the definition of pressure (7) and the fast convergence properties of the cluster expansion of P(C), we find: 

^SRB /g- E,s.^ l°g I dot LoSi\+C J2,^,^ log I det DSeoheoSi \ \ 

P(C) = lim 4t log -^—^ 5 , (9) 

where: 

1. the matrix L = C o was introduced in section 2.2 above; 

2. /U^J is the unperturbed SRB measure: if 0{tp) is a local holder continuous observable, it is defined as 

and, independently of the boundary conditions, it is equal to the Lebesgue measure. 
Defining Uq as 

Uc = log I det i| - C log I det S^^ oDSeohe\, (11) 
and using that is the Lebesgue measure on TJv, we find: 

P(C) = lim — log ^ ^ — , (12) 

so that P'(0) is equal to 

P'(0) = lim ^ V /cl^log|det.,-^Z..a^e(^^^))|e-^.^-.-°^-^-^^ _ 



1 , E 



e 



-ExcA(^i-c<^^)(-^) 
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Since P'(0)|f=o if' trivially = 0, we can try to see if d^P' {0)\e=Q is different from zero (if it were, P'(0) would be 
different from zero for £ ^ small enough). Recalling that f{ip) is the perturbing function and f^{ip) is its projection 
on the ^-th site, we get: 



Since / is periodic we have d^P' {0)\e=o = 0. 
A straightforward calculation shows that 



~Tr[{So'Df{Sli;)f] - ^ TT{S^'DfiSli;))TT(-\S^'Df{S^,i;))} , (15) 



j'eiT 



where Tr ^"^ is the trace restricted to the (unperturbed) unstable manifold. The preceding expression can be rewritten 
in a more convenient way. Using the explicit expression of (5/i(i) eq.(5), and defining Aq = Uj£„„(o)nn(^), we find that 
the first term in eq.(15) is equal to 

EE/ 7^%^A-/°"'°"«"^(V')(-a2)A-+--/€°=(C*'^'"''"^V) = 

i«i<i ~ 

EE/ (^^[l^A-r^'"" Wa2A-/^"^'«'"^(C^"^'"''"^V) . (16) 



»i = ± p>0 

iei<i ~ 



Integrating by parts, we see that the sum of the second and third term in eq. (15) is equal to: 

Combining the three contributions, we finally find: 

d',P\0)U^o = - E E/ 7^t^A-/°"'°°\^)A-/«°^'^"^(5o^^) = 



(18) 



that is the expected Green-Kubo relation (see [GR97]). 

/^From cq.(18) and eq.(8), we see that, for e small enough, 77+ is negative and, generically, strictly negative (the 
condition for / to be generic is just that the first line in eq.(18) is different from 0). 

Let us now compute eq.(18) in a special case, essentially the simpler possible. Let 

f\^)= E M^l-^v) ' f'W=0- (19) 

?7Gnn(J) 
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Substituting such choice in eq.(18), we find: 



where ei = (1, 0) and we used that v+ = {^j^^, —^j^^)- 

So, choosing e e R small enough and different from zero, r]+ = P'{0) 
S M has modulus smaller than 1, P(C) is strictly convex (since ^P"{0) = 
and £ small enough). 

B.3 Lcirge deviations 

In the present section we shall prove a large deviations property for ijaq ■ We will follow the classical strategy set up 

in [Si], [GLM02] (in particular wc will refer to the formulas in section 5 of the latter). The proof below will hold in 
the case ?7+ < 0, namely in the generic case or, to be definite, in the case the perturbation is chosen as in eq.(19). 

Thanks to the convexity of P{C), given r] € [P'(— 1), P'(l)], there exist a unique point Z{r]) G [—1,1] such that 
P'(^Z{rj)) = rj; considering such a point i] and its neighbor of radius S, Is{ri)'> such that Isiv) C [P'(— 1), P'(l)], from 
the "large deviation property III", section 5 of [GLM02], we get: 

/x"^^(j7+ e 75(77)) =0(l)e°(*l^''l)eO(l^^''l)exp{[P(Z(77)) -P(0)-Z(?7)?7]|Ao|} . (21) 

In our case P(0) = 0. Still for 77 G [P'(— 1), P'(l)], we define the free energy F{r]) as the Laplace transform of the 
generating function P(C): 

F{n)=^ max{Cr? - P(C)} = Z{rj)rj - P{Z{r,)) ; (22) 
therefore, for Is{ri) C [P'(-l), P'(l)] 

/x^'^=^(r7+ e IM) = 0(l)eO(*l^°l)eO(l^^°l)exp{ - |Ao|AF(77, r?+)} ; (23) 

where AF{r],ri+f= F{r]) - F(r/+) (indeed P(?7+) = -P(0) = 0). 

Finally, if [a,b] C [P'(-l), P'(l)], it is suitable to take (5a„ = |Ao|^'^, < /3 < 1, and divide the interval [a,b] in 

\b — a||Aop identical disjoint subintervals centered in rirf= a + (n — 1/2)(5ao, ^I'^d: 

|fe-a||Ao|'' 

/x^-=^(7?Aoe[a,6]) = M'""(^+e/5Ao(^n)) (24) 

n=l 

= 0(l)|Ao|^eO(l^<'l'"')eO(l^^<'l)exp||Ao| max f - AP(77, 7?+)l | ; 

L »)e[a,6] J 

namely the result in the second theorem. 
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= -jqifrrs-e^ + 0{s^) < 0. Furthermore, if 
-P'(0) > and P(C) is analytic for |C| < 1 
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